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ABSTRACT 

We present Submillimeter Array (SMA) A = 0.88 mm observations of an in- 
frared dark cloud G28.34+0.06. Located in the quiescent southern part of the 
G28.34 cloud, the region of interest is a massive (> 10^ M©) molecular clump PI 
with a luminosity of ~ 10'^ L©, where our previous SMA observations at 1.3 mm 
have revealed a string of five dust cores of 22-64 Mq along the 1 pc IR-dark fila- 
ment. The cores are well aligned at a position angle of 48° and regularly spaced 
at an average projected separation of 0.16 pc. The new high- resolution, high- 
sensitivity 0.88 mm image further resolves the five cores into ten compact con- 
densations of 1.4-10.6 Mq, with sizes a few thousands AU. The spatial structure 
at clump (~ 1 pc) and core (~ 0.1 pc) scales indicates a hierarchical fragmen- 
tation. While the clump fragmentation is consistent with a cylindrical collapse, 
the observed fragment masses are much larger than the expected thermal Jeans 
masses. All the cores are driving CO (3-2) outflows up to 38 km s~^, majority 
of which are bipolar, jet-like outflows. The moderate luminosity of the PI clump 
sets a limit on the mass of protostars of 3-7 Mq. Because of the large reservoir of 
dense molecular gas in the immediate medium and ongoing accretion as evident 
by the jet-like outflows, we speculate that PI will grow and eventually form a 
massive star cluster. This study provides a first glimpse of massive, clustered 
star formation that currently undergoes through an intermediate-mass stage. 

Subject headings: ISM: individual (G28.34+0.06) — ISM: jets and outflows — 
stars: formation — stars: early-type 
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1. Introduction 

Stars arc born in dense molecular clouds through gravitational collapse. A classic picture 
of low-mass star formation has been described in Shu et al. (1987). During the collapse, a 
parsec-scale molecular clump fragments into gaseous cores that subsequently collapse and 
fragment to condensations which eventually form individual or a group of stars. After the 
formation of an accretion disk, a bipolar wind develops to shed angular momentum in the 
disk so that circumstellar matter can accrete onto the protostar. Interactions between such 
a wind and the ambient molecular cloud may produce highly coUimated molecular outflows 
(or "molecular jets") at parsec scales (Shang et al. 2006, 2007). 

Massive protostars are often found embedded in massive molecular clumps of 10^ Mq 
(Plume et al. 1997; Mohnari et al. 2000; Beuther et al. 2002b; Beltran et al. 2006). Their 

presence is marked by high luminosity (> 10^ -^o), copious spectra from organic molecules 
(Blake et al. 1987; Schilke et al. 1997), active accretion (Zhang 2005; Cesaroni et al. 2006), 
and strong molecular outflows (Zhang et al. 2001, 2005; Beuther et al. 2002c). However, the 
stage prior to the high mass protostellar phase is not clear. Do massive protostars evolve 
through a low to intermediate-mass stage? How does accretion proceed? Do cores harboring 
massive young stellar objects collect all the mass before the formation of a protostar? These 
questions can potentially be addressed by studying massive infrared dark clouds (IRDCs). 

IRDCs are dense (> lO^cm"^) and cold (< 25 K) molecular clouds that absorb Galactic 
background mid-infrared emission while emit at longer wavelengths (> 100 /xm). Because 
of their extreme properties, IRDCs are ideal laboratories to study the early stages of star 
formation, and thus have been subjected to intense studies since their identification in mid- 
1990s by the Infrared Space Observatory (ISO) and the Midcourse Space Experiment {MSX; 
Perault et al. 1996; Egan et al. 1998; Hennebelle et al. 2001). Comprehensive catalogs have 
been established based on Galactic IR surveys from the MSX (Simon et al. 2006a,b) and 
Spitzer (Peretto & Fuller 2010), revealing more than 10,000 IRDCs. Among them, the ones 
with masses of 10^ Mq and luminosities < 10^ Lq are candidates of massive star formation 
ge prior to the high mass protostellar phase. 

G28.34+0.06 (hereafter G28.34), one of the flrst recognized and well studied IRDCs 
(Carey et al. 1998, 2000), appears to be a candidate to form massive stars. The giant 
cloud, at a distance of ~4.8 kpc, contains several 10^ Mq of dense gas traced by 1 mm dust 
emission along the IR-dark filaments extending ~6 pc (Carey et al. 2000; Simon et al. 2006b; 
Rathborne et al. 2006; Pillai et al. 2006; Wang et al. 2008; Zhang et al. 2009). Two prominent 
dust clumps, PI and P2, are revealed from dust continuum images obtained from single dish 
telescopes (Carey et al. 2000; Rathborne et al. 2006). Large amount of dense gas (~ 10^ Mq 
within < Ipc) makes PI and P2 potential sites of protostellar cluster formation. Zhang 
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et al. (2009) observed PI and P2 with the SMA at 230 GHz band. They found that, despite 
similar amount of dense gas contained in the two clumps, PI and P2 show very different 
stages of evolution. At a ~ 1" resolution, the northern clump P2 fragments to two compact 
cores with masses 49 and 97 Mq, respectively. The cores emit rich organic molecular line 
emissions analogous to that observed in hot cores. The southern clump PI, on the other 
hand, fragments to a string of five cores regularly spaced along the IR-dark ridge, with masses 
of 22-64 Mq. In contrast to the cores in P2, the PI cores show no molecular hne emission 
except faint C0(2-l) across the entire 4 GHz SMA band. This implies a significant CO 
depletion at a scale over 0.07 pc in the PI cores. In addition, an enhancement of deuterium 
fractionation was also found toward PI (Chen et al. 2010). All these properties point PI to 
a quiescent stage much earlier than P2. However, strong H2O mascr (Wang et al. 2006) and 
enhanced 4.5 ;um emission ("green fuzzy"; Chambers et al. 2009), as well as a faint 24 //m 
source in PI suggest that star formation may be ongoing. 

Here we report our new SMA observations of G28.34-P1 at the 345 GHz band. The high 
angular resolution, high sensitivity observations reveal hierarchical structures and multiple 
jet-like outflows, which suggest that intense star formation is already ongoing in this region. 

2. Observations and Data Reduction 

We observed the molecular clump PI in IRDC G28.34-I-0.06 with the SubmiUimeter 

Array ^ (SMA; Ho ct al. 2004) at the 345 GHz band from 2009 April 26 through 2009 
September 23. During six observing runs, the SMA antennas were in three different array 
configurations: compact (COM), extended (EXT), and subcompact (SUB), each was used 
to observe for two nights for this project. Table 1 summarizes the observations. 

For all the observations, we used a common phase center 
(R.A., dccl.) J2000 = (18M2"'50."74, -04°03'l5."3). Quasars J1751+096 and J1830+063 were 
observed approximately every 20 minutes to monitor time dependent antenna gains. The 
receivers were tuned to a local oscillate frequency of 341.6 GHz, with a uniform spectral 
resolution of 0.812 MHz (or 0.7km s"^) across the entire band. The system temperature 
varies from 120 to 500 K during the six tracks. The full-width-half-maximum (FWHM) 
primary beam is about 34". 



^The SubmiUimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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For the COM observations, we used quasar 3C273 for frequency dependent bandpass 
calibration. Absolute flux was obtained by comparing the observed correlator counts with 
modeled fluxes of dwarf planet Ceres. With an IF of 4-6 GHz, the lower sideband (LSB) 
covered rest frequencies from 335.7 through 337.7 GHz and the upper sideband (USB) covered 
345.6-347.6 GHz. 

For the EXT observations, 3C454.3 was used for bandpass calibration, while young star 
MWC349 or planet Uranus was observed for flux calibration. The observations made use 
of the newly upgraded capability that doubled the original IF bandwidth to 4 GHz in each 
sideband, covering rest frequencies 333.7-337.7 GHz in the LSB and 345.6-349.6 GHz in the 
USB. 

For the SUB observations, bandpass was calibrated by observing 3C454.3 and flux was 
scaled by observing Jovial satellites Callisto and Ganymede. The frequency coverage was 
the same as the EXT observations. 

The visibility data were calibrated using MIR IDL software package.^ Calibrated visi- 
bility data were then exported to MIRIAD^ for further processing and imaging. CASA^ was 
also used for part of the image processing. 

Data from different observing runs were calibrated separately, and then combined in 
the visibility domain for analysis. Continuum emission was generated by averaging line free 
channels in the visibility domain. In this study, if not otherwise stated, we use the EXT 
data for continuum analysis and the combined SUB and COM data for spectral analysis. 
With natural weighting, the synthesized beam is 0".69 x 0".64, with PA = — 83°.l for the 
continuum image; while for the spectral line, the beam is 3". 23 x 1".96, PA = —47°. 4. The 
la rms noise is about 0.8 mJy in the continuum, and about 100 mJy (0.2 K) per 1km in 
the spectra. If we use the SUB data only, the beam becomes larger and we obtain a la rms 
of 0.1 K at the same spectral resolution. Empirically, the absolute flux is accurate to about 
15%. The uncertainty in the absolute position is < 0".l, obtained by comparing the fitted 
position of the secondary gain calibrator and its catalog position. 

^ http:/ /www. cfa.harvard.edu/^cqi/mircook. html 
^ http:/ /www. cfa.harvard.edu/sma/miriad 
^ http://casa.mao.edu 
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3. Results and Discussion 

3.1. Continuum data and fragmentation 

3.1.1. Hierarchical fragmentation 

There are different definitions of chimp, core, and condensation in the hterature when 
describing the spatial structure of dust continuum or molecular line emission. In the following 
discussion, we refer a clump as a structure with a size of ~ 1 pc, a core as a structure with a 
size of ~ 0.1 pc, and a condensation as a substructure of ~ 0.01 pc within a core. A clump 
is capable of forming a cluster of stars, a core may form one or a small group of stars, and a 
condensation can typically form a single star or a multiple-star system. This nomenclature 
is consistent with the one adopted by Zhang et al. (2009). 

Figure 1 presents the 0.88 mm continuum emission. The 0".69 x 0".64 resolution image 
resolves the clump into five groups of compact condensations, corresponding to the five dust 
cores discovered at 1.3 mm, namely SMAl, SMA2, SMA3, SMA4, and SMA5 (Zhang et al. 
2009). Cores SMAl, SMA3, and SMA5 consist of one condensation. Cores SMA2 and 
SMA4 consist of three condensations, with one dominating over the other two in peak flux. 
In addition, a faint condensation is also revealed northeast to core SMA2. We tentatively 
assign this condensation to the SMA2 group. We name all the 10 condensations accordingly 
and label them on Fig. 1. Refining the core positions using the positions of the dominant 
condensations, the five cores are spaced by an projected distance of 0.16±0.02pc (6".8±1".0), 
and are well aligned at a position angle (PA, east of north) of 48°. In cores SMA2 and SMA4, 
the three condensations are spaced by 0.03 ± 0.007 pc (1".3 ± 0".3). 

We estimate the mass of each condensation. Assuming optically thin dust emission, the 
dust mass can be estimated following 



-"'-'dust /rp \ ) 

^ u[ J- dust )l^u 

where Afjust is the dust mass, F^, is the continuum flux at frequency u, d is the source distance, 
-Biy(Tdust) is the Planck function at dust temperature Tdust, and = 10(z//1.2 THz)^ cm^ g^^ 
is the dust opacity (Hildebrand 1983). Wc adopt the dust temperatures from Zhang ct al. 
(2009), i.e., 16 K for SMA2 and 13 K for other cores. The dust opacity at millimctcr/sub- 
millimeter wavelengths is uncertain. For dust in the interstellar medium, /3 = 2 (Draine & 
Lee 1984). Rathborne et al. (2010) derived /3 — 1.5 ± 0.3 for the entire PI clump based on 
a global spectral energy distribution. We generated a dust opacity map by comparing the 
continuum image at 345 GHz and 230 GHz, and found /3 to vary from 1.3 to 2.5 across the 
map. Given the uncertainties, we adopt (3 — 1.5 for all the cores similar to Zhang et al. 
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(2009), and discuss the effect of different ^ on the mass estimates. 

Table 2 hsts coordinates, size, peak flux, integrated flux, and estimated mass of all the 
10 condensations. These parameters are determined by fitting a two dimensional Gaussian 
function to the observed flux distribution. Fluxes are then corrected for primary beam 
attenuation. The masses of the condensations range from 1.4 to 10.6 with an average of 
4.8 The mean size is 0".9 x 0".7, or 0".6 x 0".2 after deconvolved with the synthesized 
beam. The deconvolved size corresponds to about 3000 x 1000 AU at the source distance. 
The uncertainty in the mass estimation arises from several factors. First, the flux calibration 
contributes 15 % in uncertainty. Second, the uncertainty of /3 contributes in a form of 
Mdust oc S.S'^. Third, the spatial filtering effect of radio interferomctry preferentially picks up 
the compact structures and filters out the extended emission. Comparing images made from 
different array configurations, the EXT image recovered 20 %-30 % fiuxes of what covered by 
COM and SUB images. This filtering effect explains why the masses reported in Table 2 are 
consistently lower than those in Zhang et al. (2009). The EXT image at 345 GHz represents 
the high contrast structure in these cores — the small, compact condensations tracing the 
immediate surroundings of the protostars. 

Our observations reveal fragmentation at different spatial scales in G28.34-P1. First, 
the 1 pc clump fragments into five cores (Zhang et al. 2009); second, two of the 0.1 pc cores 
fragment into even small condensations of 0.01 pc. The fragmentation at clump and core 
scales are consistent with a hierarchical fragmentation picture. The clump fragmentation is 
likely the result of initial physical conditions (density, temperature, turbulence, and mag- 
netic fields), while the core fragmentation results from an increased density after the initial 
fragmentation. Based on the Very Large Array (VLA) NH3 (1,1) and (2,2) spectra (Wang 
et al. 2008) and the IRAM 30m 1.3 mm continuum image (Rathborne et al. 2006), we mea- 
sure an initial gas temperature of 13 K and density of 3 x 10^ cm~^, over a scale of 1 pc 
toward the PI clump. These values yield a thermal Jeans mass of 0.8 M©. At the core scale, 
density increases to several 10 times of the initial clump density. For instance, in SMA2 
the core density measured from the 230 GHz image is 9.6 x 10^ cm^'^, 32 times of the initial 
clump density. This density, together with an increased gas temperature (16 K), leads to a 
thermal Jeans mass of 0.2 Mq. These parameters are also listed in Table 3. 

It is worth stressing that, in both clump and core fragmentation, the observed fragment 
(core/condensation) masses are significantly (~ 10 times or more) larger than the expected 
thermal Jeans masses (Table 3), indicating that thermal pressure is not dominant in the 
fragmentation processes. Instead, turbulence may play an important role in supporting these 
large core/condensation masses, as discussed in Zhang et al. (2009) (see further discussion 
in Section 3.1.2). The la rms sensitivity in Fig. 1 is 0.8 mJy, corresponding to 0.2 M©, well 
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below the thermal Jeans mass at the clump scale (0.8 M©). However, we did not detect 
any cores of 0.8 Mq in the filament. The non-detection of Jeans mass cores in the clump 
may suggest that low mass cores are not significantly centrally peaked, thus they cannot be 
substantiated from the smooth emission in the clump. 

3.1.2. Cylindrical collapse 

The most intriguing feature in G28.34-P1 is the configuration of well aligned, regularly 
spaced cores. The alignment and regularity strongly suggest a mechanism that is responsible 
for shaping the molecular filament into what it is now. Similar features (but at different 
spatial scales) have been reported in a large number of nearby dark clouds (Schneider & 
Elmcgrccn 1979), in a few filamentary IRDCs (e.g., Jackson ct al. 2010; Micttincn & Harju 
2010), as well as in numerical simulations (Martel et al. 2006). It has been suggested that 
these fragments are most hkely the results of gravitational collapse of a cyhnder, as first 
proposed by Chandrasekhar & Fermi (1953) and followed up by Nagasawa (1987) (also see 
discussion in Jackson et al. 2010). The theory predicts that, under self gravity, the gas in 
the cylinder will ultimately beak up into pieces with a certain interval (Amax) at which the 
instability grows the fastest and thus dominates the fragmentation process. Such 'sausage' 
instability produces a chain of equally spaced fragments along the filament, with the spacing 
roughly the interval Amax- The mass per unit length along the cylinder (or linear mass 
density) has a maximum value, (M/Omax = 2vVG = 465 {^^) MqPc-\ where G is the 
gravitational constant. If the cylinder is supported by thermal pressure, v is the sound speed 
Cg] if, on the other hand, it is mainly supported by turbulent pressure, then v is replaced by 
the velocity dispersion a. 

We now compare the observational results with the theoretical predictions. In G28.34- 
Pl, the FWHM fine width is 1.7 km s"-*^, measured form the NH3 (1,1) spectrum (Wang et al. 
2008; Zhang et al. 2009). The fine width is equivalent to a velocity dispersion of cr = 0.72 
km s~^ if the line profile is Gaussian. This velocity dispersion leads to a maximum mass 
density of 240 Mq pc~^. The total core mass in PI is 183 Mq over 0.8 pc along the filament 
(Zhang et al. 2009), implying a linear mass density of 230 Mq pc~^. Taking into account the 
missing fiux in the interferometer maps, the observations match the theoretical predictions 
well if G28.34-P1 is mainly supported by turbulence. 

The separation between two adjacent fragments, Amax, depends on the nature of the 
cylinder. For an incompressible fiuid cylinder in absence of magnetic field, Amax = 11^, 
where R is the unperturbed cylindrical radius. And the instability manifests at about two 
free-fall timescales (see Section 3.2.2). In G28.34-P1, the average separation between cores 
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is Amax = 0.16 pc, implying an initial radius oi R — 0.014 pc, several times smaller than the 
core diameter of 0.1 pc (Zhang et al. 2009). If magnetic field is involved, the strength is 
of the order of -B = AwpRy/G (see Equation (93) in Chandrasekhar & Fermi 1953). With 
n = 3 X 10^ cm~^ and R = 0.014 pc, we obtain a magnetic field strength of 0.16 mG. This 
is relatively small compared to magnetic field strength observed in massive star formation 
regions (e.g., Girart et al. 2009). 

For an isothermal gas cylinder, Amax = 22i/, where H = v{AnGpc)^^^'^ is the scale 
height, whereas Pc is the gas density at the center of the cylinder. The equation can be 
simplified as 



In the above equations, the sound speed Cg of 0.21 km s is calculated using a gas temper- 
ature of 13 K; volume density 3 x 10^ cm~^ is the averaged clump density, while 5 x 10^ 
cm~^ is the average of clump density (3 x 10^ cm~^) and core density (9.6 x 10^ cm~^) 
given in Table 3. We see that both scenarios are in agreement with the observed spacing of 
Amax = 0.16 ± 0.02 pc if the assumptions are reasonable. However, since the density should 
be the central density of the cylinder, we suspect that the turbulent support is more likely 
dominant in G28.34-P1. 

Comparison between observational results and theoretical predictions suggests that the 
fragmentation in the G28.34-P1 filament is well represented by a cylindrical collapse. The 
filament is likely supported mainly by turbulence rather than thermal pressure. 



We measure the core density profiles in unresolved cores SMAl, SMA3, and SMA5. The 
measurements are performed in the visibility domain to avoid defects in image deconvolution. 
We combine the SUB, COM, and EXT data for analysis for this purpose. The combined 
visibility data sample the G28.34-P1 region over baseline lengths ranging from about 5 to 
200 kA, offering a fairly well sampled {u, v) measurement. 

For a spherical core, the observed fiux distribution over radial distance F{r) represents 
the radial density profile modified with temperature gradient. The dust temperature scales 
as Tdust ^ (Scoville & Kwan 1976) with a = 0.33, if the core is internally heated. As 
long as Rayleigh- Jeans approximation holds, this temperature gradient contributes to the 
total fiux in the form of oc r~". Assuming a power law density profile, p oc r~^, and optically 




for thermal support, 
for turbulent support. 



3.1.3. Core structures 
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thin dust emission, the observed dust continuum flux integrated along the hne of sight, 
F oc / pTdustds, or 

F oc r-^^+''-^') , 

when (a + b) > 1. In the visibihty domain, this is Fourier transformed into a form of 

A ry q(«+^-3) 

where Auv is the visibihty amphtude and Suv — y^Jv^~+v^ is the (u, v) distance (Looney 
et al. 2000; Zhang et al. 2009). 

Fig. 2 plots the amplitude versus {u,v) distance for SMAl, SMA3, and SMA5. The 
amplitude is a vector average of the visibility data in a concentric annuli defined by a bin 
of 8kA. Least-squares fitting yields b = 2.09 ± 0.09 for SMAl, 1.97 ± 0.08 for SMA3, and 
1.84 ± 0.10 for SMA5. The results are consistent with Zhang et al. (2009) measured from 
the 230 GHz continuum data, where they derived b = 2.1 ± 0.2 for SMA4. The density 
profiles are similar to that of an isothermal sphere in hydrostatic equilibrium, where the 
radial density scales as oc r~^. We note that the large error bars and discrepancy at long 
{u,v) distances reflect sensitivity limits and may be also partly due to unresolved weak 
sources in the vicinity. 



3.2. Spectral line data and CO outflows 

3.2.1. Morphology 

Among the entire 8 GHz SMA band, only GO (3-2) is detected above 3a {la 0.1 K at 1 
km resolution). Nevertheless, GO (3-2) reveals high- velocity outflowing gas in PI. Fig. 3 
presents the GO (3-2) channel maps, and Fig. 4 plots the integrated blueshifted emission 
in blue contours, and redshifted emission in red contours, superposed on the continuum 
emission. One can find in Fig. 3 that the extended GO emission close to the systemic 
velocity (78.4 km s~^) is filtered out. This effect actually helped us identify outfiows at 
high velocities. It is intriguing that all the five dust cores are associated with GO outflows, 
majority of which are bipolar, jet-like outflows. 

The outfiows are centered on the dust cores and arc generally oriented cross the major 
axis of the filament. The most prominent outfiow is centered on core SMA2, with a blue 
lobe emanating toward southeast (SE; see channels of 51-71 km s~^ in Fig. 3), and a red 
lobe toward northwest (NW; 86-116 km s~^). In addition, a minor blue lobe is also seen at 
NW (66-76 km s~^). The outflow centered on SMA3 is in an east-west orientation, with its 
blue lobe emanating westbound (51-71 km s~^) and its red lobe eastbound (86-96 km s"^). 
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SMA4 lies in the center of a NW-SE outflow, with a SE blue lobe (66-76 km s"^) and a 
NW red lobe (86-91 km s~^). SMAl is associated with a quadrupolar outflow. Its two blue 
lobes are seen at 66-71 km s~^, emanating toward SE and north, respectively. Its two red 
lobes arc seen from 86 through 91 km s~^ toward NW and south, respectively. The outflow 
centered on SMA5 has an orientation almost parallel to the SMA4 outflow. This outflow 
is weak compared to others, and its SE blue lobe and NW red lobe can be barely seen at 
66-71 km and 86-91 km respectively. Orientation of all the outflows are sketched 
as arrows in Fig. 4. 

The driving sources of these outflows can be idcntiflcd geometrically. Cores SMAl 
to SMA5 are roughly located at the geometric centers of the outflows, respectively. While 
SMAl, SMA3 and SMA5 are unresolved, SMA2 and SMA4 are resolved into three condensa- 
tions. In these two cases, we assign the strongest and/or central condensation as the driving 
source of the relevant outflow. Hence, the outflows are likely driven by condensations SMAl, 
SMA2a, SMA3, SMA4a, and SMA5. 

The CO outflows shown in Figs. 3 and 4 consist of pairs of knots which are geometrically 
symmetric with respect to the central protostar. Illustrated in Fig. 4, outflows SMA2a and 
SMA3 have four pairs of knots, and outflow SMA4a has two pairs. The knots are equally 

spaced by about 0.16 pc. These knots in outflows may represent outbursts due to disk 
variability arising from episodic, unsteady accretion (e.g., HH 211, Lee et al. 2010; HH 80- 
81, Qiu & Zhang 2009). With the spacing and maximum outflow velocity fmaxj we estimate 
the period between two outbursts to be about 4000 yr. This value lies in the range of roughly 
10^ — 10'^ yr used in episodic accretion models for low mass stars (Baraffe & Chabrier 2010; 
Zhu et al. 2010). 

The outflows show slight bending (Fig. 4). In the outflow SMA3, the red lobe bends 
5° toward the south, and the blue lobe bends 6° toward the south, making the outflow C- 
shaped. In outflows SMA2a and SMA4a, and probably SMAl, the knots seem to trace the 
propagation of a periodically wiggled, S-shaped molecular jet, most evident at the closest 
knotty pairs near the protostar. Part of this effect may be contaminated by multiple outflows 
(e.g., in outflow SMA2a). Nevertheless, the mirror symmetric (C-shaped) wiggles and point 
symmetric (S-shaped) wiggles may indicate different bending mechanisms like jet precession, 
orbital motion of the jet source, and Lorentz forces (Fendt & Zinnecker 1998; Masciadri & 
Raga 2002; Raga et al. 2009). Similar bending effects have been observed in a number of 
Herbig-Haro jets (e.g., HH 211, Lee et al. 2010) and some of the most collimated high-mass 
molecular jets (e.g., Su et al. 2007). Because the wiggles are not prominent in CO (3-2), we 
do not discuss this effect further in this paper. 
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3.2.2. Physical properties 

We calculate the physical parameters of each outflow, including mass, momentum, en- 
ergy, as well as dynamical age and outflow rate. Assuming local thermodynamic equilibrium 
and optically thin CO emission in the line wings, we first derive the CO column density 
following Garden et al. (1991), 



N, 



/33 12\ f 

^co(cm-2) _ 4 ^ ^Qi2 (y^^ ^ Q l-j^j J TBdv , 

where dv is the velocity interval in km s~^, and Tex and Tb are excitation temperature and 
brightness temperature in K, respectively. We take the gas temperatures to be the excitation 
temperatures, i.e., 16 K for the SMA2a outflow and 13 K for other outflows. The outflow 
mass, momentum, energy, dynamical age, and outflow rate are then given by 

M = d^ [^]mn, J^Ncoin')dn' , 

P^Mv, 

E = -Mv"^ , 
2 

, _ -^flow 
^dyn ) 

11 

'^max 

Mont = :: — , 

^dyn 

where d is the source distance, fn-n^ is the mean mass per hydrogen molecule assumed to be 
2.33 atomic units, Q is the total solid angle that the flow subtends, v is the outflow velocity 
with respect to the systemic velocity of 78.4km s~^, Vmax is the maximum outflow velocity, 
and Lflow is the flow length. We adopt an empirical [^] abundance ratio of lO'' (Blake et al. 
1987). 

Table 4 lists derived physical parameters for each outflow. These outflows show sym- 
metry in blue- and redshifted lobes. The energetics are small compared to outflows observed 
in high mass protostellar objects (Beuther et al. 2002c; Zhang et al. 2005), but are similar 
to outflows emanating from nearby (< 1 kpc) intermediate-mass hot cores (IMHCs: IRAS 
22198+6336, Sanchez-Monge et al. 2010; NGC 7129-FIRS 2, Fuente et al. 2005; IC 1396 
N, Neri et al. 2007). Particularly, the outflow rate Mout, which amounts to (4 — 47) x 
lO~^M0yr~^, and the flow momentum P, which ranges from 1 to 9 AfQkm s^^, are of the 
same order to all the three IMHC outflows known to date. The dynamic timescales are 
(1.5 — 3.4) X 10"^ yr for all the outflows. These values are comparable to massive outflows but 
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are one order of magnitude higher than the IMHC outflows. Because observations are hmited 
by noise and gas at higher velocities are possible, tdyn is a lower limit of the outflow age and 
the accretion history of the dominant protostar. Theoretical models suggest that the massive 
star formation process takes a few times of the free-fall timescale {ts— ■^■^^ — ^J7==|=) 

of the core. Beuther et al. (2002c) observed a rough equality between tdyn and tg. For the 
PI cores, we obtain tg of (1.2 — 6.7) x 10'^ yr based on the densities given in Table 3. We see 
that tfj approximates tdyn well, and this supports the idea that flow ages are good estimates 
of protostar lifetimes (Beuther et al. 2002c). 

The SMA2a, SMA3, and SMA4a outflows are so well coUimated that they resemble 
"molecular jets" commonly found in low-mass star-forming regions (e.g., L1157, Zhang et al. 
2000; OMC-IS, Zapata et al. 2005; HH 211, Lee et al. 2007). Recently, interferometric 
observations have revealed collimated outflows in high-mass star-forming regions (e.g., AFGL 
5142, Hunter et al. 1999; Zhang et al. 2007; IRAS 20126+4104, Cesaroni et al. 1999; Shepherd 
et al. 2000; IRAS 05358+3543, 19410+2336, 19217+1651, Beuther et al. 2002a, 2003, 2004; 
HH80-81, Qiu & Zhang 2009). To measure the coUimation, we take the ratio of outflow total 
length (blue + red lobes) over width measured in Fig. 4 as coUimation factor. We obtain 14, 
10, and 7 for SMA2a, SMA3, and SMA4a outflows, over scales of 0.5-1 pc. These values are 
comparable to that of the most collimated low-mass outflows and are higher than that of 
the high-mass outflows discovered by single dishes and previous interferometers (Wu et al. 
2004; Beuther et al. 2002c; Qiu & Zhang 2009). Beuther et al. (2007) reported one outflow in 
IRAS 05358+3543 with a coUimation factor of ~ 10 and claimed that it was the first massive 
outflow observed with such a high degree of coUimation at a scale of 1 pc. In G28.34-P1, the 
SMA2a and SMA3 outflows are among the best coUimated molecular outflows yet discovered 
in high-mass star-forming regions over parsec scales. We note that these coUimation factors 
should be considered as lower limits because the outflow lengths and widths are not corrected 
for inclination angle, and are not deconvolved with synthesized beam. 



3.3. Massive star formation through an intermediate-mass stage 

The presence of highly collimated outflows demonstrates a protostellar nature of the 
condensations in G28.34-P1. Wang et al. (2008) reported a luminosity of 10^ for the 

24 fim source associated with SMA2. Rathbornc et al. (2010) obtained a luminosity of 
2.1 X 10'^ Lq for the entire PI clump. These luminosities set a limit on the mass of the 
embedded stars of 3-7 Mq, if they are at the zero-age main sequence. The stellar masses 
put these objects at an intermediate-mass stage. The presence of highly collimated outflows 
also favors a disk-mediated accretion scenario of massive mass star formation (Zhang et al. 
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2005; Cesaroni et al. 2007) in a manner similar to low- mass star formation. In addition, the 
high coUimation also indicates that these outflows are at very early stages since the opening 
angle of outflows widens as they evolve. 

The G28.34-P1 condensations are comparable in mass to the nearby intermediate-mass 
"cores" (2-5 Mq) at a similar spatial extent (2300 x 1600 AU; Sanchez- Monge et al. 2010). 
The non-detection of organic molecules (CH3OH, CH3CN, etc.) in PI is likely due to the 
sensitivity limit of our observations. Sanchez- Monge et al. (2010) observed copious hot core 
molecular line emissions toward 122198, an intermediate-mass star with a luminosity of 370 
Lq at a distance of 0.76 kpc. The typical brightness temperature of the lines is 0.8 K. If 
we scale 122198 to the distance of G28.34, the emission would become 40 times fainter, 
well below our detection hmit of 0.1 K per 1 km s~^ resolution (also true for Zhang et al. 
2009). Therefore, it is possible that G28.34-P1 may have already experienced an increase 
in abundance of organic molecules as expected from protostellar heating, similar to 122198. 
Deep spectroscopy with ALMA will be able to verify the presence of hot core type molecules 
in G28.34-P1. 

Despite the similarities between the G28.34-P1 condensations and the intermediate- 
mass star forming "cores" , there is an important distinction between the two: The G28.34- 
Pl condensations arc embedded in compact (~ 0.1 pc), dense cores (~ lO^cm"^, sec Table 
3) with several tens of solar masses, which are hierarchically embedded in a massive, dense 
(3 X 10^ cm~^. Table 3) molecular chimp of over 10^ M©, an order of magnitude larger than 
the mass in a typical intcrmcdiatc-mass star-forming region. With a typical star formation 
efficiency of 30% in dense gas and a standard stellar initial mass function, it is natural to 
expected that massive stars will eventually emerge in some of the condensations (if not all) , 
at the end of protostellar accretion. At that time, PI would become a massive star cluster. 
In the context of evolution, an immediate comparison to G28.34-P1 is the northern clump 
G28.34-P2 in the same cloud complex, which may represent the future of PI. With a similar 
clump mass, P2 has a higher gas temperature of 30 K, and a higher luminosity of 2.6 x 10"^ 
Lq (Rathborne et al. 2010), equivalent to a zero-age main-sequence mass of 14 Mq. The 
SMA observations detected molecular lines from CH3OH and CH3CN in the P2 cores, while 
none in the PI cores (Zhang et al. 2009). Likewise, other massive molecular clumps (e.g., 
W33A, Galvan-Madrid et al. 2010; G31.41, Girart et al. 2009; NGC 6334, Hunter et al. 2006; 
G10.6, Liu et al. 2010) that have luminosities of 10^ — 10^ Lq, and are associated with radio 
continuum emission, may represent an even more evolved phase of massive star formation 
at which high mass protostellar objects have already emerged. 

The large mass reservoir (22-64 at ~ 0.1 pc; 1000 M© at < Ipc) in G28.34-P1 and 
the outflows suggest that the condensations are, and will remain, actively accreting materials 
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through the ambient medium. Materials in the cores fall onto the embedded protostars 
through disks and at the same time materials in the clump fall onto the cores. Observations 
of G28.34-P1 seem to indicate a continuous accretion over a spatial scale that spreads two 
orders of magnitude, from 1 to 0.01 pc, and finally reaches onto the protostars. This is a 
relatively new picture, but has also been reported in another IRDC, G30.88+0.13 (Zhang 
& Wang 2011). A comparison with massive and more luminous molecular clumps indicates 
that the condensations will grow in mass, and likely emerge as massive stars eventually. The 
embedded protostars appear to be still intermediate-mass objects and are on their way to 
eventual massive stars. Our observations of the G28.34-P1 clump provide perhaps the first 
glimpse of an intermediate-mass stage of massive, clustered star formation. 



4. Conclusion 

We present an observational study of the PI clump in IRDC G28. 34-1-0. 06 using the 
SMA at the 345 GHz band. Our main findings are summarized as follows. 

1. High resolution and high sensitivity image at 0.88 mm reveals hierarchical fragmen- 
tation at Ipc and 0.1 pc scales. The mass of the dust cores and condensations are much 
larger than the relevant thermal Jean masses. The clump fragmentation is consistent with 
a cylindrical collapse. 

2. All five dust cores are driving highly coUimated CO (3-2) outflows, indicating a 
protostellar nature of the condensations embedded in the cores. Bipolar, jet-like outflows 
suggest a disk-mediated accretion of high-mass star formation. 

3. The relatively low luminosity and a lack of detection of organic molecules indicate 
that the condensations in PI underline an intermediate- mass stage of massive star formation 
and will likely form massive stars eventually. 
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Tabic 1. Summary of the SMA Observations 



Array 


UT date 


Antenna 


Bandwidth 


Bandpass 


Flux 


T 


T225GHZ ^ 




(yyyymm) 


No. 


(GHz) 


Calibrator 


Calibrator 


(K) 




COM 


20090426 


7 


2 


3C273 


Ceres 


200-400 


0.09-0.15 


COM 


20090503 


7 


2 


3C273 


Ceres 


160-280 


0.05 


EXT 


20090806 


8 


4 


3C454.3 


MWC349 


260-500 


0.09-0.14 


EXT 


20090826 


7 


4 


3C454.3 


Uranus 


120-250 


0.06-0.07 


SUB 


20090905 


5 


4 


3C454.3 


Ganymede, Callisto 


160-300 


0.05-0.07 


sue 


2()()9()92.S 


t 


4 


.'-',r4.l4.:-! 


Gaiuiiu'dc. Callisto 


i.')0-:-!00 


0.07-0.09 



'^Zenith opacity measured from water vapor monitors mounted on CSO or JCMT. 



Table 2. Properties of the Condensations 



Name 


RA(J2000) 


Dec(J2000) 


Size 


P.A. 


-fpeak 


-Pint 


Mass/M© 




li:m:s 


d:m:s 


a" 


X b" 


O 


mJy 


mJy 


(3 = 1.5 


/3 = 2 


SMAl 


18:42:51.19 


-04:03:07.2 


0.82 


0.68 


-63.6 


25.3 


31.6 


10.6 


19.8 


SMA2a 


18:42:50.85 


-04:03:11.4 


0.74 


0.70 


55.2 


15.7 


18.3 


4.3 


8.1 


SMA2b 


18:42:50.76 


-04:03:11.5 


1.18 


0.68 


56.8 


8.5 


15.5 


3.6 


6.8 


SMA2c 


18:42:50.79 


-04:03:12.4 


0.81 


0.62 


60.5 


5.1 


5.8 


1.4 


2.6 


SMA2d 


18:42:50.99 


-04:03:09.5 


0.94 


0.63 


62.2 


5.1 


6.7 


1.6 


3.0 


SMA3 


18:42:50.58 


-04:03:16.3 


0.76 


0.63 


-82.3 


14.4 


15.6 


5.2 


9.8 


SMA4a 


18:42:50.28 


-04:03:20.2 


0.78 


0.75 


38.1 


12.2 


16.0 


5.3 


10.0 


SMA4b 


18:42:50.32 


-04:03:21.0 


1.30 


0.75 


33.2 


6.8 


15.0 


5.0 


9.4 


SMA4c 


18:42:50.20 


-04:03:20.3 


0.93 


0.55 


-84.7 


7.4 


8.5 


2.8 


5.3 


SMA5 


18:42:49.83 


-04:03:25.2 


0.91 


0.74 


78.8 


16.6 


25.1 


8.4 


15.7 
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Table 3. Clump/ Core Fragmentation 



Source 


rj-i a 


n(H2) 


Mj 


Rj 


Fragment Mass ^ 




K 




Me 


pc 


Mo 


PI 


13 


3.0 X 10^ 


0.8 


0.04 


22-64 


P1-SMA2 


16 


9.6 X 10^ 


0.2 


0.008 


> 1.4-4.3 



^Gas temperatures are derived from VLA NH3 (1,1) and 
(2,2) lines (Wang et al. 2008). 

^Observed core/condensation masses are calculated using a 
dust opacity index of /3 = 1.5. 



Table 4. Derived Outflow Parameters 



Parameter 


SMAl 


SMA2a 


SMA3 


SMA4a 


SMA5 




Blue 


Red 


Blue 


Red 


Blue 


Red 


Blue 


Red 


Blue 


Red 


V (km s~^) 


[67,76] 


[85,93] 


[41,76] 


[84,116] 


[50,71] 


[85,99] 


[64,72] 


[85,107] 


[66,73] 


[86,93] 


M (Mq) 


0.20 


0.07 


0.38 


0.23 


0.22 


0.18 


0.06 


0.12 


0.05 


0.06 


P (Makm s-i) 


1.12 


0.71 


5.09 


4.04 


3.17 


2.05 


0.52 


1.50 


0.35 


0.63 


E (Mq km2 s-2) 


3.85 


3.62 


53.46 


46.18 


26.96 


13.06 


2.57 


11.15 


1.40 


3.40 


iflow (pc) 


0.40 


0.27 


0.51 


0.48 


0.56 


0.51 


0.44 


0.43 


0.34 


0.39 


tdyn (10^ yr) 


3.43 


1.81 


1.33 


1.25 


1.93 


2.42 


2.99 


1.47 


2.68 


2.61 


Mout(10-5MQyr-i) 


0.58 


0.41 


2.86 


1.84 


1.14 


0.74 


0.20 


0.82 


0.19 


0.23 



''Parameters are not corrected for inclination angle. 
''Range of velocities used to derive outflow parameters. 
'^Projected length of the longer lobe in case of two lobes. 



J2000 Right Ascension 



Fig. 1. — SMA 0.88 mm image of the G28.34-P1 region in contour overlaid on the Spitzer 
24 fim image in color scale. The SMA image is made from the EXT configuration data only, 
with synthesized beam 0".69 x 0".64, PA = — 83°.l, indicated as an ellipse in the lower left 
corner of the panel. The contours start at 4mJy (Sex) and increase by a step of 1.6 mJy {2a). 
Assigned condensation names are also labeled on the image. The SMA image shown here is 
not corrected for the primary beam attenuation. 
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Fig. 2. — Amplitude distribution against {u,v) distance of cores SMAl, SMA3, and SMA5. 
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— 111.0 
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(7 
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50 [- 106.0 
0300 
lO' - 
20" 
30" 
-04°03'40" - 



— 116.0 




O °* SMAl 
: OJ *SMA2a^ 

O (^MA4a o- 
*:SMA5 - 



RA (J2000) 

Fig. 3. — CO (3-2) channel maps. The image is made by averaging the emission over 
a 5 km s~^ velocity interval, with the central velocity labeled on each panel. The contours 
start at 5(7 (240 mJy) and step by ±5(T. The filled ellipse on the bottom left corner represents 
the synthesized beam. The stars mark the presumable outflow driving sources as labeled on 
the last panel (see the text). 




^gh42m52^5 51 ^5 51 ^0 50^.5 50^0 49^5 49^0 

J2000 Right Ascension 

Fig. 4. — Integrated CO (3-2) emission overlaid on the 0.88 mm continuum emission of 
the G28.34-P1 region. The velocity ranges used for integration are [41, 74] km s~^ for the 
blueshifted lobe (blue contours), and [84, 116] km s~^ for the redshifted lobe (red contours). 
The CO contours are ±(2, 2.5, 3, 9.5) x 10% of the blue lobe peak, 39 Jy beam~^-km s~^. 
The continuum contours start at 4mJy (5cr) and increase by a step of 1.6 mJy {2a). The 
white cross represents the H2O maser detected by VLA (Wang et al. 2006). [Follow up 
observations have resolved this maser spot into two, spatially coincident with SMA2a and 
SMA2b, respectively (Wang et al. 2011 in prep.)] The arrows sketch outflow directions. The 
synthesized beams indicated in the lower left corner are 3". 23 x 1".96, PA = —47°. 4 for the 
CO image, while 0".69 x 0".64, PA = — 83°.l for the continuum image. 
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